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ABSTRACT 



We compute the mass and composition of dust produced by stars with masses 
in the range IMq ^ M ^ SMq and with a metalhcity oi Z = 0.001 during their 
Q ■ AGB and Super AGB phases. SteUar evolution is followed from the pre-main sequence 

' phase using the code ATON which provides, at each timestep, the thermodynamics 

C/2 , and the chemical stucture of the wind. We use a simple model to describe the growth 

of the dust grains under the hypothesis of a time-independent, spherically symmetric 
stellar wind. Although part of the modelling which describes the stellar outflow is not 
completely realistic, this approach allows a straight comparison with results based on 
^ ' similar assumptions present in the literature, and thus can be used as an indication of 

CO \ the uncertainties affecting the theoretical investigations focused on the dust formation 

l/^ . process in the surroundings of AGB stars. 

' We find that the total mass of dust injected by AGB stars in the interstellar 

, medium does not increase monotonically with stellar mass and ranges between a min- 

imum of 10~^Mq for the l.SM© stellar model, up to 2 x 1O~^M0, for the 6M0 case. 
Dust composition depends on the stellar mass: low-mass stars (M < SM©) produce 
carbon-rich dust, whereas more massive stars, experiencing Hot Bottom Burning, 
never reach the carbon-star stage, and produce silicates and iron. This is in partial 
disagreement with previous investigations in the literature, which are based on syn- 
thetic AGB models and predict that, when the initial metallicity is Z = 0.001, C-rich 
' dust is formed at all stellar masses. The differences are due to the different modelling of 

H \ turbulent convection in the super-adiabaticity regime. Also in this case, like for other 

physical features of the AGB, the treatment of super-adiabatic convection shows up 
as the most relevant issue affecting the dust-formation process. 

We also investigate Super AGB stars with masses O.SM© ^ M ^ %Mq that evolve 
over a ONe core. Due to a favourable combination of mass loss and Hot Bottom 
Burning, these stars are predicted to be the most efficient silicate-dust producers, 
releasing [2 — 7] x 1O~^M0 masses of dust. 

We discuss the robustness of these predictions and their relevance for the nature 
and evolution of dust at early cosmic times. 

Key words: Stars: abundances - Stars: AGB and post-AGB. ISM: abundances, dust 



1 INTRODUCTION 

Winds from stars on the asymptotic giant branch (AGB) and 
super-asymptotic giant branch (SAGB) provide an impor- 
tant component of mass return into the interstellar medium 
(ISM) and account for a significant fraction of interstel- 
lar dust in present-day mature galaxies (Zhukovska, Gail 
& Trieloff 2008; Sloan et al. 2009). Recent chemical evo- 



lutionary models, which attempt to constrain the origin of 
the large dust masses inferred from mm and submm data 
of 2; > 6 QSOs, have shown that AGB stars give important 
contributions to the dust content even at these early cosmic 
epochs (Vahante et al. 2009, 2011). 

Observations with the Spttzer telescope have probed 
the formation of dust in AGB stars of galaxies in the Lo- 
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cal Group (Sloan et al. 2008, 2010a and references therein) 
and in globular clusters (Sloan et al. 2010b; McDonald et 
al. 2011). In spite of continuous progress, detailed under- 
standing of the outflow dynamics and grain formation is far 
from being complete (for a comprehensive review see Hofner 
2009): dust grain properties depend on dynamical aspects of 
the stellar wind, and on the surface chemical composition of 
the star. Shock waves caused by stellar pulsations lift the gas 
above the stellar surface intermittently creating dense cool 
layers where grain condensation can occurr. Provided that 
their cross sections are large enough, the newly formed dust 
grains can be accelerated by radiation pressure and collision- 
ally couple to the gas, dragging it along. Time-independent 
wind models, though limited in their capability of reproduc- 
ing important properties of the dynamical atmosphere with 
shocks, can reproduce certain observed features, such as the 
observed spectral energy distribution of C-rich stars with 
high mass loss rates. 

Compared to the problem of computing reliable mass 
loss rates, the problem of calculating dust yields has the ad- 
ditional complication of estimating the fraction of elements 
condensing into various types of dust species, which, in turn, 
depends on the changes in the surface chemistry of the star 
determined by the interplay between nucleosynthesis and 
convection. 

The most extensive study so far has been made by Fer- 
rarotti & Gail (2006). Using synthetic stellar evolution mod- 
els, they compute the non-equilibrium dust formation and 
estimate dust yields for M-, S-, and C-type AGBs. Syn- 
thetic AGB models are useful tools, due to the relatively 
short computation time required and the simplicity with 
which they can be incorporated in stellar population syn- 
thesis studies. On the other hand, they have no predictive 
power on some key physical processes that occurr during 
AGB evolution, such as the Hot Bottom Burning experi- 
enced at the bottom of the convective envelope. A proper 
description of these processes and of their impact on dust 
formation demands the full integration of the stellar struc- 
ture. 

To progress in this direction, the aim of the present in- 
vestigation is to estimate the mass and composition of dust 
formed around AGBs using models that follow stellar evo- 
lution from the pre-main sequence phase. This enables the 
computation of the time evolution of the surface chemistry 
and of the mass loss rate in a self-consistent way, provid- 
ing an important benchmark for studies based on synthetic 
modelling. Using the proper stellar evolutionary sequences, 
we also apply the dust formation model to winds of SAGB 
stars, whose dust yields have not been studied so far. 

The paper is organized as follows. The description of the 
codes used to calculate the AGB and SAGB evolution, the 
wind structure and dust formation are described in Sect. 2; 
the main physical processes which affect the surface chem- 
istry of these stars are discussed in Sect. 3. In Sect. 4 we 
present the resulting dust mass and composition for stars 
with masses in the range IMq ^ M ^ SM© and with initial 
metallicity of Z — 0.001; we assess the robustness of these 
results and discuss the effects of uncertain physical processes 
in Sect. 5. Finally, in Sect. 6 we compare our predictions with 
alternative dust yields present in the literature and in Sect. 
7 we summarize the main conclusions. 



2 PHYSICAL AND NUMERICAL INPUTS 

To address the issue of dust production by AGB stars, we de- 
veloped a specific tool to calculate the structure of the wind 
around the central object as well as the mass and compo- 
sition of the newly formed dust. These properties depend 
on the luminosity, effective temperature and mass loss rate 
experienced by the star, and on the surface chemical com- 
position. In the following, we provide a description of the 
numerical code used to integrate the stellar evolution, the 
ATON code, and of the code used to estimate the stellar 
wind structure and dust formation. 



2.1 The ATON evolution code 

The evolution of the stars, from their pre Main Sequence 
phase, until the almost complete ejection of their exter- 
nal mantle, wa s computed by means of the code ATON 
(|Ventura et al.l il998) which integrates spherically symmetric 
structures in hydrostatic equilibrium. An exaustive descrip- 
tion of the code is beyond the scope of the present study and 
we refer the interested reader to the original paper. Here we 
briefly recall those features which affect the thermodynamics 
and chemical composition of the wind where dust nucleation 



2.1.1 The convection model 

The temperature gradient within regions unstable to con- 
vection is determined via th e Full Spectrum of Turbulence 
ICanuto fc Mazzitelli (l9_9l|), hereinafter FST) model for 
turbulent convection; the code allows the alternativ e Mix- 
ing Length Theory (MLT) treatment (| Vitensd IT953I ') . Mix- 
ing of chemicals within nuclearly active regions is handled 
by means of a diffusive approach: for each of the species in- 
cluded in the nuclear n etwork we solve the diffusive equation 
following the scheme bv lCloutman fc Eolll l|l97(j ). which cou- 
ples self-consistently nuclear burning and mixing of chem- 
icals. The extra-mixing (overshooting) from the convective 
borders is simulated by an exponential decay of velocities 
from the neutrality point, where buoyancy vanishes; the 
scale of this phenomenon is assumed to be / = ^Hp, where 
Hp is the pressure scale height. In agreement with the cal- 
ibration based on the fit of the observ e d ma in sequences 
of open clusters given in IVentura et al] (|l998l l. we assume 
C = 0.02 to simulate overshoot from the border of the con- 
vective cores during the two major phases of core hydrogen 
and helium burning. To limit the number of free parameters, 
we assumed no extra-mixing from the base of the convective 
envelope during the giant evolution, and during the AGB 
phase. This choice renders the extent of the Third Dredge 
Up found during the Thermal Pulses phase a conservative 
estimate. 



2.1.2 The mass loss rate 

The mass loss rate du ring t he AGB evolution was modelled 
according to jBlock3 jl99i). Based on hydrodynamic com- 
putations by Bowen (|l988l ). the strong increase in the mass 
loss rate during the AGB evolution is modelled by multiply- 
ing the Reimers' rate by a given power of the luminosity; 
the resulting expression is. 
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2.2 Stellar wind and dust formation 



M = 4.83 X IQ-^^tirM-^'^U-^R 



(1) 



where riu is a free parameter and M, R, and L are the stellar 
mass, radius and luminosity, expressed in solar units. Using 
as a calibration the luminosity function of lith ium-rich stars 
in th e Magellanic Clouds, we set rjn = 0.02 (jVentura et al.l 
l2000l ). 

The Blocker treatment has been used and tested so far 
to describe mass loss from high-luminosity, M stars; here we 
extend this modelling to the C-rich regime, being aware of 
the additi onal uncertainti e s due to this choice on the results 
obtained. iMattsson et al.l (|2007l ) present an analysis of the 
uncertainties in the description of the mass loss during the 
AGB phase of carbon stars (see in part icular their Fig. 4). 

Some sequences based on the IVassiliadis fc WoodI 

l|l993l) treatment (hereinafter VW93) were also calculated 
for comparison: in this case mass loss increases exponen- 
tially with the pulsation period, P, until the star enters a 
super-wind phase, when P exceeds 500d. The super-wind 
rate is M = 5 x lO'^Me/yr. 

The results provided by the two descriptions are in rea- 
sonable agreement for low-mass AGBs, and for the early 
AGB phases of the more massive objects; the main differ- 
ences are found in the description of the advanced phases 
of the evolution of the sta rs with mass c lose to the limit for 
carbon ignition, where the lBlockeij (|l995l ) treatment predicts 
very large mass loss rates, that favour an earlier consump- 
tion of the whole stellar envelope. 



The thermodynamic and chemical structure of the wind that 
forms around AGBs is calculated on the basis of the inputs 
provided at each time step by the integration of the struc- 
ture of the central star, which are the natural outcomes of 
the ATON code. In particular, we use luminosity (L), effec- 
tive temperature (Teff), mass (M), mass loss rate (M), and 
surface chemistry defined by the mass fractions (Xi) of the 
various elements included in the nuclear network. 



2. 2. 1 The structure of the wind 

The structure of the wind is det e rmine d following the 
schematization by iFerrarotti fc Gaill (|2006l ). 

The outflow is assumed to be stationary and spherically 
symmetric. Neglecting the hydrostatic pressure forces, we 
find from the equation for momentum conservation that. 



dv 
dr 



GM 



(l-P). 



(2) 



The quantity P gives the ratio between the radiative pressure 
on the dust and the gravitational pull and is given by. 



kL 



where k is the flux-averaged extinction coefficient of the 
gas-dust mixture and can be expressed as. 



y ] fi 



(4) 



2.1.3 Nuclear network and opacities 

The nuclear network includes 30 elements, with 64 re- 
actions. All the most relevant p-capture and a-capture 
processes are taken into account. The nuclear cross-sections 
were taken from the NACRE compilation (jAngulo et al.l 
Il999l). with the excep tion of the rea ctions ^'^N(p 7)^''0 
■ '^Cra.'Y) "^0 (iKunz et al.[ |2002| ). 
C+7 (|Fvnbo et all [iooi ), ^^Ne(p , 7)^^N a 



3a 

l|Hale et al.l |2002|) 

^^Na(p.a)^°Ne l|Hale et a"Lll2004h 



Na (p,7r ''Mg (|Hale et al 
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J mi), 



The opacities in the low-temperature regime are cal- 
culated by means of th e AESOPUS tool, developed by 
iMarigo fc Aringed l|2009l ). These tables are suitably con- 
structed to follow the changes in the chemical composition 
of the envelope when the C/O ratio approaches, or exceeds, 
unity. This is deeply different from the classic approach, tra- 
ditionally used to build AGB models, that neglects possible 
variations in the surface chemistry during the evolution in 
the opacity calculations. To understand how the results de- 
pend on this choice, we decided to calculate additional AGB 
models, where the above effect is neglected. It goes without 
saying that this is expeted to play a role only in the models 
that eventually reach the C-star stage, i.e. for M< 4Mq in 
the present investigation. 

For the present study, we restrict our discussion to 
stellar models with absolute metallicity Z=0.001, and 
helium fraction Y=0.24. The mixture is a-enhanced 
with [a/Fe]=-|-0.4 , the solar reference mixture being 
iGrevesse fc Sauval l|l998l ). Other chemical compositions will 
be investigated in future explorations. 



with itgas = IQ-^p'^^^T^ l|BeU fc Linlll994l ). The sum in 
eq. (|4]) is extended to all the dust species considered: the fi 
terms give the degrees of condensation of the key-elements 
for each dust species and ki represent their corresponding 
extinction coefficients. 

To completely define the thermal structure of the wind, 
we need to specify the density and temperature radial pro- 
files. Mass conservation leads to, 



M — 4Tvr pv, 



(5) 



and the temperatu re structure i s dete rmined by adopting the 
approximation bv lLucvl (|l97ll , ll97d) that holds for spheri- 
cally symmetric, grey winds: 



1 -2^cff 



1- A 1- 



R2 



(6) 



The optical depth r in eq. Q is found from the differential 
relation, 



dr 
dr 



— —pk 



R^ 



(7) 



with the limiting condition that r — > for r — > 00. 

The wind model, determined via integration of eqs. ([2|- 
(O, is calculated starting from the innermost radius where 
the first dust species become stable, commonly located 3-4 
stellar radii away from the centre of the star. This is ex- 
pected to coincide with the region where P increases until 
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it exceeds unity, and the wind is accelerated to s upersonic 
velocities. In analogy with lFerrarotti fc Gaill (|2006l ). we keep 
the velocity of the gas constant from the surface of the star 
to the region where dust forms. The above equations are 
integrated out to a radial distance of 10^*7? {R is the stellar 
radius), far beyond the point where all the relevant quanti- 
ties reach their asymptotic behaviour. 

Before turning to the dust formation process, we stress 
here that a more physically sound description of the wind 
properties is needed before the full reliability of the results 
from this kind of investigations can be confirmed. In particu- 
lar, the assumption of a time-independent flow neglects any 
effect of atmospheric shocks waves, which play a crucial role 
for the wind mechanism. The assumption of a constant ve- 
locity does not account for the consequences of these shocks, 
which could push the gas beyond the condensation point. 



thermal velocity, reading 



^ OliVth. 



Pv,i 



(10) 



For carbon, based on the discussion in iFerrarotti fc Gail 
(2006), we assume that condensation begins at a temper- 
ature of 1100 K and that the destruction term is negligible. 

The vapour pressures pv,i in Eq. 10 are found via the law 
of mass action: it allows to express the individual equilib- 
rium gas pressures of the molecules involved on the basis of 
the total variation of enthalpy associated to the dust forma- 
tion reaction. An exhaustive description of the methodology 
followed to es timate the individual vap our pressures can be 
found, e.g., in IFerrarotti fc Gaij l|200ll ). 

The formation of the four dust species in M-stars (C/O 
< 1) is therefore characterized by the following four equa- 
tions, 



2.2.2 Dust formation 

Dust formation generally starts with the condensation of 
small seed nuclei which then grow by accreting material to 
macroscopic grains. Since the nature of these seeds is still 
unclear, we do not consider the nucleation phase, but only 
the growt h process of grains . As in itial conditions, in agree- 
ment with lFerrarotti fc GaiJ (|2006l ). we assume that the seed 
grains have an initial size ao = Inm, and that their initial 
density is n^j = 3 x 10~''nH, where uh is the hydrogen num- 
ber density in the wind. We could verify that the results are 
independent of both assumptions. 

The composition of dust in the winds during the AGB 
phcise depends on the surface chemical composition. For M 
stars, with surface C/O < 1, the dominant dust species are 
expected to be olivine, pyroxene, quartz, and solid iron. For 
C stars with surface C/O > 1, the dust mixture is assumed 
to be dominated by solid carbon, SiC and solid iron. In Ta- 
ble 1 we summarize the formation reactions of these dust 
species. When computing the reaction rates for olivine and 
pyroxene, we follow the variation of the magnesium percent- 
age with respect to the pure iron component. 

The temporal variation of the dust grain size for the 
i-th dust species, a^, is computed as a competition between 
the growth rate Jf' and the destruction rate Jf^'^, i.e. 



(8) 



where Vo,i is the volume of the nominal molecule in the 
solid. For each dust species, the growth rate is determined 
by the addition rate of a key element, which is generally the 
least abundant species involved in the chemical reaction. For 
silicates and SiC these are supposed to be the SiO and C2II2 
molecules, respectively, whereas for solid iron and carbon it 
is the corresponding element in the gas phase. The growth 
rates per unit time and surface area are proportional to the 
thermal velocity, Vth, and density, n, of the key species in 
the gas phase, 



(9) 



where the index i runs over the different key species and a 
is the sticking coefficient. 

Vaporization of dust grains by thermal decomposition 
is estimated from the vapour pressure of the key species, 
Pv,i, over the solid state, and it is also proportional to the 



dt 



Vo,qu( J; 



(11) 



(12) 



(13) 



(14) 



These equations are completed by two further relations, giv- 
ing the percentages of magnesium in the olivine and pyrox- 
ene dust particles: 



dt 
dt 



3Vo,o 



{Xg — Xpy) Jpy + — (J+ — J_ ) 



(15) 



(16) 



where Xg is defined as the relative gas abundance of magne- 
sium with respect to the Mg+Fe sum. 

The quantity {J^ — Jl") gives the difference between 
the exchange rate of iron by magnesium per unit surface 
area during collisions of Mg with the grain surface and the 
rate of the reverse reaction, and is given by 



(J™ — J°^) = Vth,MgCt^'^ i^riMg ~ UirKp) 



(17) 



where Kp is calculated on the basis of the free enthalpies 
of formation of the pure products with only magnesium or 
iron, and q*^^ is the exchange coefficient, describing the prob- 
ability of exchange of magnesium and iron in the scattering 
process. 

Dust formation in C-stars is described by eq. (|14p . plus 
two additional equations for carbon and SiC grains. 



C^Q-SiC T7 / TQ-r 7dcc\ 



dt 



dac ,gr 



(18) 



(19) 
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Table 1. Dust species considered in the present analysis, their formation reaction and adopted sticking coefficient (see text). 
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Grain Species Formation Reaction Sticking Coefficient 



Olivine 2a;Mg +2(l-a;)Fe+SiO+3H20 ^- Mg23:Fe2(i_3;)Si04 + 3H2 0.1 

Pyroxene xMg +(l-a;)Fe+SiO+2H20 -J- Mga;Fe(i_^)Si03 + 2H2 0.1 

Quartz SiO + H2O -> Si02(s) +H2 0.1 

Silicon Carbide 2Si + C2H2 2 SiC + H2 1 

Carbon C ^ C(s) 1 

Iron Fc — > Fc(s) 1 



The growth rate coefficients J^l, J|y, J|Jj, J?', 



and the exchange coefficients J^, Jl" for olyvine 
and p yrox ene are calcii l a,ted fo llowing Gail k, Sedlmavij 
d 1999V Tcrra rotti fc Gaill (l200lD . iFerrarotti fc Gaill (|2002l ) 
and iFerrarotti fc Gaill (|20oi r The sticking coefficients for 
iron, carbon and SiC formation are assu med to be q = 1, 
where as for olivine we take a = 0.1 (see IFerrarotti fc Gaill 
l|2002l ). Table 3). For pyroxene and quartz the sticking co- 
efficients are less well determined from l abora tory experi- 
ments; in analogy with IFerrarotti fc Gaill (|200ll ). we assume 
Qpy = Qfqu = 0.1 (see the last column of Table 1). Finally, 
the exchange coefficient entering the equation for the deter- 
mination of the fraction of magnesium within olivine and 
pyroxene is taken to be oi'^ = 0.06. 

The free enthalpy of formation of the molecules involved 
in the reactions for the condensation of the vario us species 
of dust were taken from ISharp fc Huebneil (|l990l ). with the 
only exceptions of FeSiOa, graphite, and solid SiC, that were 
kindly provided by Prof. Gail (private communication). 

The extinction coefficients for iron and quartz were 
calculated according to the analytic expansion given in 
iGail fc Sedlmavil (| 19991 ). whereas for olivine and p y roxen e 
we used the opacities derived bv lOssenkopf et all l|l992h . 
For what concerns the C-rich env ironment, we use for SiC 
eq. (34) in IFerrarotti fc Gaill (|2002l ). whereas the carbon ex- 
tinction coefficient was kindly provided by Prof. Gail (pri- 
vate communication). In the cases where analytic expan- 
sions associated to solar abundaces were used, the coeffi- 
cients were properly scaled, to account for the difference in 
the number density of the key-species involved in the present 
investigation, compared to the solar values. 



3 THE CHANGE IN THE SURFACE 

CHEMISTRY OF AGB AND SAGB STARS 

After the core He-burning phase is finished, stars with mass 
1M(7 ) ^ M < Mu p evolve through the AGB phase (see 
e.g. iHerwid (|2000l )): they experience a series of thermal 
pulses (TP), triggered by the periodic ignition of a He- 
rich layer below the CNO burning shell, and end-up as CO 
White Dwarfs. More massive objects, with Mup <M<MccSn, 
achieve carbon burning in an off-center, partially degenerate, 
region; after experiencing some TPs, they evolve as White 
Dwarfs made up of oxygen and neon. Their evolution is com- 
monly referred to as the Super Asymptotic Giant Branch 
(SAGB) phase. 

The exact values of Mup and MccSn are determined by 
the core overshooting during the two major phases of nuclear 



burning; in the present investigation, based on the assump- 
tion of a moderate extra-mixing from the outer border of 
the convective core, we find Mup = 6M0 and MccSn = 8M0. 

The mass and composition of dust produced by AGBs 
and SAGBs is determined both by the physical behaviour 
of the star and by the surface chemistry. The amount of 
dust produced will depend on the number density of the key 
species, which demands a detailed understanding of the evo- 
lution of the surface mass fraction of silicon, oxygen, carbon, 
magnesium and iron. 

The surface chemistry of AGBs is affected by two fun- 
damental physical process es, the Third Dr edge Up and the 
Hot Bottom Burning (e.g. Ilben fc Renzinil ([l983)). In what 
follows, we briefly describe the effects of these two mecha- 
nisms, and their associated uncertainties. 



3.1 Third Dredge Up 

The Third Dredge Up (hereafter TDU) is associated to the 
inwards penetration of the surface convective zone following 
each thermal pulse (TP), when the bottom of the external 
mantle crosses the entropy barrier due to the H-He discon- 
tinuity; the envelope enters a region where a pulse-driven 
convective zone, formed as a consequence of the 3a reactions 
ignition, had spread the products of helium nucleosynthesis; 
the outcome of such a mechanism is surface carbon enrich- 
ment, and the possible achievement of the C-star regime. 
The TDU modelling is extremely sensitive to the numer- 
ics adopted (|Straniero et al.lll997l ): a narrow temporal and 
spatial zoning, with time-steps of the order of hours, and 
~ 5000 mesh-points, are required to allow a self-consistent 
description of the TDU phenomenology. A deeper TDU can 
be obtained either by assuming some o vershoot from the 
bottom of the surface convective zone l|Herwig|[2000l ). or 
from the borders of the convective shell that f orms when 
the thermal pulse starts (|Herwig fc Austin|[20o3 ). 

For C-rich models, a further source of uncertainty is 
the adopted low-T radiative opacities. When the C/0 ratio 
approaches unity, the main absorbers of radiation switch 
from the oxygen bearing molecules, such as H2O, VO, TiO, 
to C-bearing molecules, such as CN, C2, C3, C2H2, which 
are much more sensitive to the surface mass fractions of 
the main elements: this makes the C-rich models much 
more sensitive to the details of the variation of the surface 
chemical composition of the star. The increase in the 
absorption of radiation favoured by a larger carbon content 
is commonly neglected in the computations of the AGB 
phase, that assume the opacity to depend only on the 
initial metal content, and disregard any modification of the 
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Figure 1. The evolution of tiie luminosity (Left), mass loss rate (Middle, expressed in M0/yr) and effective temperature (Right) with 
the total mass of the star for 4 models with initial mass 6M0 (dashed— dotted), 4Mq (dotted), SMq (solid), and 2M0 (dashed). The low 
surface temperatures reached by the lower mass models are favoured by the general expansion of the whole convective envelope, that 
follows the achievement of the C-star stage. 
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Figure 2. The evolution of the surface chemistry of models with different initial masses (same symbols as in Fig. [TJ as a function of the 
stellar mass. The three panels refer to the surface abundances of Carbon (Left), Oxygen (Middle) and total Magnesium (Right). The 
arrows indicate the initial abundances of the same elements in the original mixture. We note the transition from Llot Bottom Burning, 
operating in the high-mass regime, and determining a depletion of all the three elements, to the effects of the III Dredge-up, outlined 
by the increase in the surface carbon mass fraction, dominating in the lowest masses. 



surface chemistry. As initially suggested bv iMarigol (|2002| ). 
and later con f irmed on the basis of full AGB comput ations 
l|CristaUo al.l l2008l : IVentura fc Marigol |2009| . l20ld ). the 
usage of the correct opacities favours a sudden expansion of 
the structure once C/O reaches unity, which is accompanied 
by an increase in the mass loss rate experienced; the 
consumption of the envelope is consequently faster, which 
prevents the possibility that enormous amounts of carbon 
are accumulated in the surface regions. An interesting 
consequence, to be explored also within the context of the 
present investigation, is that the general cooling of the 
structure that follows the approach of the C-star stage may 
potentially quench the Hot Bottom Burning, and favour a 
C-rich environment. 



3.2 Hot Bottom Burning 

The description of the Hot Bottom Burning (HBB) phe- 
nomenology requires a full integration of the whole stellar 
structure, and a diffusive description of the coupling of nu- 
clear burning and mix i ng of chemicals wit hin the envelope 
(Sack mann et al.ll 19741 : 1 Ventura et al.ll2000l ). The HBB phe- 
nomenon is a consequence of the delicate coupling between 
the outer region of the degenerate core, the CNO burn- 
ing layer, and the innermost regions of the surface convec- 
tive zone, and by definition cannot be described within the 
framework of a synthetic modelling of the AGB phase. 

Unlike the TDU, HBB conditions are reached during 
the interpulse phase. The temperature at the base of the 
external mantle may exceed ~ 40 x lO^K, with the conse- 
quent activation of p-capture nucleosynthesis. The products 
of such nuclear activity, due to the rapidity of the convec- 
tive currents, are immediately transported to the surface. 
A mild activation of HBB leads to the depletion of carbon, 
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Figure 3. Left: Grain size evolution as a function of tiie total mass of the star during the AGB evolution of a model with initial mass 
6Mq. The different lines correspond to olivine (dashed, black), pyroxene (dotted, red), quartz (solid, green) and iron (dot-dashed, blue). 
Right: The corresponding evolution of the fraction of silicon condensed into silicate— type dust (solid, black), and of the iron condensed 
(dot-dashed, blue) 
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Figure 4. Left: The evolution as a function of the total mass of the star during the AGB evolution of the grain size (left) and the 
fraction of the key elements condensed (right) in a model of initial mass 4Mq . The meaning of the different lines in the two panels is 
the same as in Fig. \3\ 



whereas a stronger HBB may potentially lead to the destruc- 
tion of the surface oxygen (T> 70 x 10''K); in both cases, a 
great enhancement of the surface nitrogen is achieved. The 
modelling of the HBB phenomenology is extremely uncer- 
tain, because the thermodynamic stratification of the inter- 
nal region of the surface convective zone is critically depen- 
dent on the convective model adopt ed: a higher efficienc y 
convective description favours HBB jRenzini 
[Blocker fc Schoenberneil (|l99ll ) showed that HBB ignition 
determines a steep increase in the luminosity of the star, and 



a significant deviation from th e classic rel a tionsh ip between 
core mass and luminosity by IPaczvns^ (|l97Ci ). The ma- 
jor impact of the HBB phenomenology on the AGB evolu- 
tion was investigated bv lVentura fc D'Antonal (|2005al ): when 
convection is modelled following the FST scheme, such as in 
the reference version of the ATON code, HBB is found in all 
models for stars more massive than ^ SMq, and is accom- 
panied by a rapid increase in the mass loss rate experienced 
by the star, and thus by a faster consumption of the whole 
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envelope. Only a few thermal pulses are experienced, which 
prevents the possibility that the C-star stage is reached. 



3.3 SAGB stars 

The main physical properties of massive stars with 
Mup <M<MccSN, which ex perience the SAGB phase, have 
been describ ed in detail by ISies"3 (|2006l ). The recent full ex- 
ploration bv ISiesd (|2010l ) provides yields that allow to esti- 
mate how these stars pollute their surroundings, and to infer 
the extent of the nucleosynthesis experienced at the bottom 
of their surface convective zon e. The recent investigation 
bv lVentura fc D'Antonal (|201l|} 



. vi a a de tailed comparison 

with the results presented in lsiesi (|2010l ). showed that: (i) 
the core masses of this class of objects is so high that HBB 
conditions are easily achieved, independently of the convec- 
tion modelling; (ii) the strength of the thermal pulses ex- 
perienced is modest, thus TDU never occurs, and the C- 
star stage is never achieved; [Hi) the mass loss treatment 
plays a major role in determini ng the physical evolu t ion o f 
SAGBs. When a treatment a la I VassiUadis fc WoodI (|l993l ') 
is adopted, the star experiences a large number of TPs; the 
extent of the nucleosynthesis at the bottom of the co nvec- 
tive mantle is so high that oxygen is severel y deple ted l|Siessl 
I2OIQ V On the contrary, when the ; Blocker! (|l995l ) recipe is 
adopted, the star lose all the envelope before a strong nu- 
cleosynthesis may occur. 

The SAGB stars are thus expected to experience a 
strong mass loss since the beginning of their TP evolu- 
tion, and to show-up traces of p-capture nucleosynthesis, 
although the rapid consumption of their external mantle 
prevents the achievement of extended modifications of their 
surface chemistry. 



4 WHICH KIND OF DUST FROM AGB AND 
SAGB STARS? 

The quantity of dust formed around a star depends on the 
gas density in the region where condensation occurs. Within 
the present schematization, mass conservation imposes a di- 
rect relationship between M and p (see Eq. 5): higher densi- 
ties correspond to large mass loss rates, which renders AGB 
and SAGB stars particularly suitable for the dust forma- 
tion process in their surroundings. In the following we de- 
scribe the results of the model for stars with initial masses 
in the range IMq ^ M ^ 8M0 and initial metallicity of 
Z = 0.001. The main evolutionary properties of these models 
wer e presented and discu ssed in Ventura & Marigo (2010) 
and I Ventura fc D Antonal (|201ll ). Fig. [T] shows the evolution 
of the main physical properties of four models differing in 
their initial mass; the three panels refer to the whole AGB 
phase, during which the stellar envelope is gradually con- 
sumed, with a decrease in the stellar mass. The variation of 
the surface chemistry in the same models is shown in Fig. [2] 
We limit this description to carbon, oxygen and magnesium, 
because iron is unchanged, and silicon shows a modest vari- 
ation (~ -1-0.05 dex), and only in the 6M0model. 



4.1 Dust formed under HBB conditions 

The left panels of Figs.Oand [3]shows the evolution of dust 
grain sizes formed in winds of stars with masses M = QMq 
and M = AMq (right panel), that experience HBB. The 
right panels show the fraction of the key elements condensed 
into dust: silicon for olivine, pyroxene and quartz, and iron 
for solid iron. The three panels of Fig. [T] show the evolu- 
tion of their main physical quantities during the AGB phase 
(dashed-dotted and dotted lines, for 6M0 and 4M0 , respec- 
tively), and demonstrate the large luminosities and mass loss 
rates reached, in particular by the M — 6A1q model. 

The depletion of the surface carbon (indicated by the 
initial drop in the carbon mass fraction in the left panel of 
Fig-El prevents the C-star stage to be reached, thus allowing 
only the formation of silicates and solid iron. Instead of the 
time coordinate, which increases from left to right, in both 
Figs. [3] and 3] we show the results as a function of the mass 
of the star: this helps a better understanding of the amount 
of dust produced, because the latter depends on the mass 
lost by the star when the size of the grains is the highest. 

Given our choice of convective overshoot from the core 
during the H-burning phase, the left panel represents the 
evolution of the largest mass model which evolves through 
the AGB phase. We can distinguish different phases during 
the evolution. Dust begins to form as soon as HBB begins, 
and the mass loss increases. Olivine is the most abundant 
species and it forms close to the stellar surface; as HBB be- 
comes stronger, the mass loss rate further increases, leading 
to an increase in the amount of olivine formed but, at the 
same time, causing a decrease of the least abundant species, 
i.e. qu artz and iron. This result, in agreement with the anal- 
ysis of iFerrarotti fc Gaill (|200ll ). is due to the rapid accel- 
eration experienced by the wind as soon as the first species 
form, which decreases the density of the key species. 

The size of the olivine and pyroxene grains decreases at 
each thermal pulse, as a consequence of the temporary drop 
in the luminosity (hence of the mass loss rate) of the star. 
Quartz and iron follow the opposite behaviour, because the 
smaller mass loss reduces the differences among the amount 
of dust formed for each species. 

We see from Fig. [3] that in the 6Mq model, when the 
total mass of the star drops below ~ 4:Mq, the depletion of 
surface oxygen (clearly indicated by the slope of the dashed- 
dotted line in the middle panel of Fig. [2]) determines a 
scarcity of water molecules, required for the formation of 
the three species of silicates considered (see Table 1). Be- 
cause the formation of a molecule of olivine requires three 
H2O molecules, olivine production stops, followed by py- 
roxene (which needs only two water molecules), and quartz 
(one water molecule). From this point on, only iron grains 
are produced. 

We thus find that massive AGB models that experience 
strong HBB conditions produce a large amount of silicate- 
type dust during the first part of their AGB evolution, and 
end-up with iron dust production once the surface abun- 
dance of oxygen is strongly reduced. A word of caution con- 
cerning this latter possibility is needed here, because the 
occurrence of an iron-driven outflow might be an artificial 
effect of the grey description adopted in the present mod- 
elling, which favours smaller condensati on distances f or the 
iron grains. Simple analytical estimates (|Hofneij|2009l ). sup- 
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ported by detailed models (|Woitkel2006l ). demonstrated that 
the wavelength dependence of the corresponding grain opac- 
ities in the near-IR leads to condensation distances well be- 
yond the reach of shock waves, which means that iron cannot 
start an outflow, unless the wind is triggered by other types 
of dust. The possibility of an outflow driven only by iron is 
thus highly uncertain. 

The 4Mq model (see Fig. |4| evolves at smaller core 
masses, and can be taken as representative of a star expe- 
riencing a milder HBB. Compared to the BM© case, due to 
the smaller M, less silicates are produced in the early AGB 
phases; on the other hand, the depletion of oxygen is less 
severe (see middle panel of Fig. [21), thus silicates continue 
to form throughout the AGB evolution. These two effects 
compensate each other, thus the two models, at the end of 
the evolution, produce approximately the same quantity of 
silicates and a total dust mass of ~ (1.2 - 2) X 1O"*M0 (see 
Table 2). The lower mass loss also leads to smaller differ- 
ences among the various types of dust formed. 



4.2 Dust from C~rich environments 

Models whose initial mass is below S.SAf© hardly reach HBB 
conditions, and evolve at smaller luminosities compared to 
their more massive counterparts (see the left panel of Fig. [T]). 
The repeated TDU episodes eventually leads to the forma- 
tion of a C-star, as indicated by the great increase in the 
surface carbon content for the models with initial mass 2M0 
and 3M0 shown in the left panel of Fig. (2] This is accompa- 
nied by a switch from silicate-type dust to C-rich particles. 
In the left panels of Figs.[5]and[6]we show the grain size evo- 
lution as a function of the mass of the star for two models 
with 3M0 and 2M0 ; the right panels show the degree of con- 
densation of the key elements (essentially, carbon and silicon 
in this case). In both cases we note an initial modest pro- 
duction of silicates, that proceeds until C/O reaches unity. 
From that point on, solid carbon and SiC are formed. The 
decrement in the dimension of carbon grains during each 
interpulse period in the SM© model is an effect of HBB, 
which is completely absent in the 2M0 model. In this latter 
case, representative of all the models that show no effects of 
HBB, a large amount of carbon grains are expected to form 
around the star. Unlike their more massive counterparts, 
low-mass AGBs progressively increase their surface carbon 
content due to the repeated TDUs, and produce more and 
more dust as they evolve, until no mass is left within the 
convective envelope. The total mass of dust formed by stars 
with M ^ 3.5M0 ranges between ~ (0.01 - 7.7) x 1O"''M0 
(see Table 2). 

4.3 Dust from SAGB stars 

SAGB stars evolve at large luminosities since the flrst TPs, 
and experience a strong mass loss rate that favours an 
early consumption of their convective envelope. Fig. Qshows 
the variation of the overall fraction of silicon condensed 
into silicate-type dust as a function of the total stellar 
mass. The different lines correspond to SAGB models with 
masses in the range 6.5M ^ Af ^ SAf© published in 
IVentura" fc D'Antonall|201ll ). For comparison, the most mas- 
sive AGB model (6M0) previously discussed is also shown. 



It is evident that the mass of silicate dust grains produced by 
SAGBs is larger than in the QMq model, for two main rea- 
sons: (i) since SAGBs evolve at large core masses, their lumi- 
nosity and mass loss are extraordinarily large since the very 
first TPs and (ii) unlike the most massive AGBs, oxygen is 
not depleted extensively, because the envelope is lost before 
the nucleosynthesis experienced at the bottom of their sur- 
face mantle can reach very advanced stages. Thus, silicates 
continue to form during the whole TPs evolution because 
water molecules are always available in their surroundings. 
The total masses of dust formed by SAGBs vary in the range 
~ (3 7) X 10"'' M0 (see Table 2). 



4.4 Stardust from AGB and SAGB stars 

In Table 2 we list the total dust masses and the masses 
of individual dust species predicted by the model for the 
investigated grid of stellar masses. Similarly to other impor- 
tant AGB phenomena, the key quantity which determines 
the composition of dust formed around the star is the core 
mass, Mcore- models with small Mcore achieve only a mod- 
est production of silicates, whereas carbon is produced in 
great quantities. In the more massive objects the opposite 
situation occurs, since HBB prevents the formation of any 
carbon-rich dust species, whereas a strong production of sil- 
icates is favoured. 

Fig. [S] shows the total mass of dust as a function of 
the initial stellar mass (solid line) and the contribution of 
individual dust species (i.e. silicates, iron, carbon, and SiC). 
The results refer to AGB and SAGB models calculated with 
the reference ATON model discussed in Sect[4l 

We note a doubled-peaked distribution, with a maxi- 
mum production of carbon by low-mass AGBs and of sili- 
cates by SAGBs; stars with 3M0 < M < 5Mq provide the 
smallest contribution to dust formation, because the pro- 
duction of carbon is inhibited by HBB, which is however too 
soft to allow a large silicates production. At stellar masses 
around ~ 3Mq (this result depends on the treatment of the 
low-T molecular opacities, see Sec. 15. ip we note the transi- 
tion from carbon-rich to silicate-rich dust, which is a pure 
effect of HBB, prevailing over the TDU effects. 

The mass of silicates formed by stars in the mass range 
[4 — 6]A^0 is approximately constant: this is a consequence 
of the strong production of silicates, achieved in the early 
AGB phases in the more massive models, counterbalanced 
by the scarcity of the water molecules available for the con- 
densation process, due to the very strong and fast depletion 
of oxygen, that characterize the latest evolutionary stages. 
The SAGB models, are, under the present assumptions, the 
most efficient silicates producers, for the reasons discussed 
in Sect. 4.3. 



5 DEPENDENCE OF DUST PRODUCTION 
ON PHYSICAL PARAMETERS 

In the previous section we have presented the results 
of the dust formation model for different stellar masses. 
Here we discuss how uncertainties in the modelling of the 
AGB/SAGB evolution and in the physics of dust formation 
may affect the robustness of these findings. In particular, 
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Figure 5. The same as Fig. [S] but for a stellar model with initial mass 3Mq. Line coding for silicates and iron (whose production is 
however negligible in the present model) is the same as in Fig. |3] Light solid, magenta lines indicate the evolution of carbon grain size 
(left panel) and of the fraction of carbon condensed (right); Long-dashed, cyan lines show the variation of the SiC grain size (left) and 
of the fraction of silicon condensed into SiC (right). We note the effects of the transition from M stars to C stars, when no silicates are 
produced. 



T , I , 1 r 




Figure 6. The same as Fig.fS] but for 

we examine the role played by the adopted low-T molec- 
ular opacities, the mass loss treatment, and the chosen set 
of sticking coefficients. A discussion of the impact of the 
adopted convective model is presented in Sect. 6. 



5.1 Low— T molecular opacities 

The models presented and discussed so far (see, in partic- 
ular. Figs. [H and [6]) are based on a modern and updated 
method to determine the molecular opacities in the low 
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stellar model with initial 2M0 (right). 

temp erature regime: the AESOPUS tool l|Marigo fc Aringe3 
[2003), 

that accounts for the changes in the chemistry of the 
envelope to calculate the absorption coefficient, k, in the 
most external regions of the star. 

In the traditional approach, the evolution of the surface 
abundances of the various species is ignored: the opacity is 
determined on the basis of the initial composition of the 
star, neglecting the increase in the surface carbon abundance 
provided by TDU. 

Since previous studies of dust production by AGBs rely 
on this latter approximation, we have quantified the differ- 
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Figure 7. The total fraction of silicon condensed into silicate 
dust grains as a function of the total stellar mass in SAGB mod- 
els with initial masses of 6.5 (dashed, blue), 7 (solid, magenta), 7.5 
(long— dashed, black), and 8 Mq (dot-dashed, green). For com- 
parison, we also show the results obtained for the most massive 
AGB model with M = GMq (dotted, red line). 




Figure 8. The total mass of dust produced by AGB & SAGB 
models of different masses (solid, black line). We also show the 
total mass of the individual components, indicated as follows: 
Silicates (dot-dashed, magenta), carbon (dashed, red), SiC (long- 
dashed, green), iron (dotted, blue). 



Figure 9. The dust produced by a SM© model in which, as in 
previous works, the increase in the surface carbon is neglected 
in the computation of the low— T opacities. The dimensions of 
the different types of dust produced are indicated with the same 
symbols as in Fig. ??. The comparison with the results shown in 
the left panel of Fig. ?? allows us to understand the role played 
by the low-T opacity description: we note that in this case HBB 
is not extinguished, thus the production of carbon— rich dust is 
limited to a narrow time— interval, before strong HBB destroys 
the surface carbon available. 



ences introduced by the new method to calculate the opaci- 
ties in terms of the total mass of dust formed. We thus cal- 
culate a further set of AGB models in which the traditional 
method was used to compute k in the low-T regime. 

As long as the C/O ratio is < 1, such as in the most 
massive AGB models which experience HBB, the evolution 
is ind ependent of the opacity treatment (jVentura fc Marigcl 
I2OI0I ). We thus restrict our analysis to lower mass stars (M^ 
3Mq), whose surface chemistry is dominated by TDU. 

Fig.|5]shows the resulting grain size evolution for a 3Mq 
star; this is to be compared with the left panel of Fig. (5] that 
refers to a model with the same mass, calculated with the 
AESOPUS tool. Both models achieve an early small pro- 
duction of silicate-type dust, that proceeds until repeated 
TDU episodes allow to reach the C-star stage. From this 
point on, the predictions diverge. When the opacities are 
computed with the AESOPUS tool, carbon-rich dust con- 
tinues to be produced during the whole AGB phase. On the 
other hand (see Fig. (Qjl , when the traditional opacities are 
used, C and SiC grains are produced only within a narrow 
range of stellar mass, after which silicates production takes 
over again. This is a consequ ence of t he qu enching of HBB, 
that was early predicted bv lMarigol (|2002l ). and later con- 
firmed bv I Ventura fc Marigcl 1 2009! ): because the core mass 
of a 3Mq star is the minimum allowing HBB conditions, the 
expansion of the envelope following the increase in the sur- 
face carbon is sufficient to turn off HBB: this is favoured by 
the rapid increase of AESOPUS opacities that accompanies 
the achievement of the C-star stage. The model shown in 
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Fig. [9] does not follow this behaviour due to the use of clas- 
sical opacities; the star eventually reaches HBB conditions 
and silicates production takes over again. 

The 3M0 model just discussed lies within the nar- 
row range of masses {8M ~ O.SM©) where the differences 
caused by the two opacity treatments are most dramatic, 
affecting both the resulting grain composition and the total 
mass of dust produced: when the traditional opacitities are 
adopted, the total mass of dust produced is 7 x 10~^Mq 
(to be compared with the corresponding value in Tab. 2, i.e. 
~ 2 X 1Q~'^Mq), most of which is silicate-type dust. Stars 
with 2.5M0 do not experience any HBB and the low- 
T opacity treatment affects only the total mass of C-rich 
dust formed, which is approximately a factor ~ 2 smaller 
when the AESOPUS tool is used. In fact, the increase in 
the surface carbon abundance is accompanied by a general 
expansion of the whole structure, and by a faster consump- 
tion of the whole envelope, resulting in a smaller surface 
carbon content. 

It is important to remark that the above conclu- 
sions dep end on the adopted mass loss rate. As shown by 
IVentura fc M arigo (2010), the physical and chemical evolu- 
tion of low-mass AGB stars depend on the interplay between 
mass loss rate and molecular opacities: a smaller M favours 
HBB and thus shifts downwards the range of masses domi- 
nated by TDU, around which the formation of carbon-type 
dust occurs. 

5.2 The mass loss treatment 

The mass loss description has important effects on dust for- 
mation around AGBs. In fact, the surface chemistry of the 
star, which determines the composition of dust species, re- 
sults from the interplay among TDU, HBB, and the con- 
sumption rate of the envelope. This is particularly impor- 
tant for massive AGBs and SAGBs, because the degree of 
the p-capture nucleosynthesis experienced at the base of 
the convective zone, hence the final elemental abundances, 
depends on the mass left in the envelope at each evolu- 
tionary stage. Therefore, a major issue is the competition 
between the rates of the relevant nuclear reactions and the 
time-scale for envelope consumption. In the reference ATON 
model, t he ma ss loss treatment follows the prescription by 
I Block er (1995'). To quantify the uncertainties introduced by 
the mass loss treatment w e have explored models where the 
IVassiliadis fc WoodI (|l993t ) prescription is used. 

The evolution of the QMq model presented in Fig. [3] 
is characterized by an early phase of strong silicates pro- 
duction, followed by a later stage, during which only solid 
iron forms. When silicates production ceases, the mass left 
on the envelope depends on the ratio between the rate of 
oxygen burning at the bottom of the external mantle and 
the rate of envelope ejection: a smaller M renders this be- 
haviour even more extreme, because it increases the fraction 
of silicates-free mass lost by the star. This is illustrated in 
Fig. 1101 The left panel shows the vari ations o f the s urface 
content of oxygen and carbon when the 'Blocker' (| 19951 ) (solid 
lines) and the V assiliadis fc W ood (1993) (dotted lines) mass 
loss treatments are used. In both models, an initial phase 
of carbon burnin g occurs , follow ed by the depletion of sur- 
face oxygen. The lBlockeJ (|l995l 'l mass loss is so strong that 
oxygen burning proceeds at the same rate as surface man- 



tle consumption. Since the lVassiliadis fc WoodI (| 19931 ) mass 
loss rate is much slower, the star experiences many more 
TPs, and eventually TDU, via an increase in the surface 
^^C, favours the formation of a carbon star. However, car- 
bon dust production is negligible because surface carbon 
enrichment is inhibited by HBB. For the QMq stellar mass 
model, the mass lost by the star when it is still 0-rich is 
AM ~ IM0, to be compared to AM ~ 3M0 in the B16cker| 
(I1995i ) model. As it is shown in the right panel of Fig. IIOI 
the production of silicates is consequently smaller and iron 
is the only dust species produced for most of the evolution. 
As a result, the lower mass loss rate leads to a reduction of 
the total dust mass produced by ~ 2 dex. 

Dust formation is progressively less affected by the mass 
loss description as the mass of the star decreases, because 
stars of smaller mass never enter the phase when only iron 
is produced, thus the production of silicates never ceases. 

As discussed in Sect. 15.11 mass loss also plays a role in 
setting the stellar mass, which marks the border between 
silicates producers and carbon -dust producer s. With our 
standard choices, based on the iBlockerl (|l995l ) treatment, 
we find this threshold m ass to be M ~ 3M0. When the 
IVassiliadis fc WoodI (|l993l ) formulation is adopted, this limit 
shifts downwards to M ~ 2.5M0, because the larger mass 
left in the envelope for a given core mass favours HBB igni- 
tion. 

For SAGBs, I Ventura fc D'Antonal (|201ll ) showed that 
a mass l oss ra te increasing with the luminosity, as in the 
iBlockeJ (|l995h recipe, prevents the star from experiencing 
a very advan ced nucleosynthesis; on the contrary, in models 
based on the Vassil iadis fc WoodI l|l99j ) formulation only a 
modest amount of oxyg en is left in the envelope for most 
of the SAGB evolution (|Siess|[2O10l ). In this latter case the 
production of silicates would be smaller and limited to early 
evolutionary phases, whereas iron would be the main dust 
species produced. A detailed investigation, aimed at under- 
standing the uncertainties affecting Stardust from SAGB 
stars, is in preparation and will be presented in a separate 
paper. 

5.3 The sticking coefficients 

For each dust species, condensation occurs below a threshold 
temperature that depends on the differences in the forma- 
tion enthalpies of the various molecules entering the reaction 
associated to the condensation process. Below the threshold 
temperature, the destruction rate rapidly declines to zero, 
thus leaving the production rate to depend uniquely on the 
growth-rate coefficients. The latter, in turn, are linearly pro- 
portional to the sticking coefficients, a, which parametrize 
the efficiency of the dust formation process. 

For quartz and pyroxene the sticking coefficients are 
poorly constrained by laboratory experiments. Similarly, 
the exchange coefficient regulating the magnesium fraction 
within olivine and pyroxene is very uncertain. Since for all 
the stellar models the surface chemical composition is ex- 
tremely rich in magnesium, the latter coefficient has a neg- 
ligible effect on dust formation and the equilibria of both 
hybrid products, olivine and pyroxene, are always shifted to 
the magnesium-rich cases. 

To estimate the uncertainties associated to the sticking 
coefficients we have artificially changed by a factor 5 those 



Stardust from AGB and Super- AGB stars 13 



Table 2. Dust mass produced by AGB and SAGB models, assuming an initial metallicity Z = 0.001. The initial stellar mass M and 
the core mass Mcorc are reported in the first and second columns. The total mass of dust, and the mass of pyroxene (Mpy), quarz 
(Afqu), solid iron (A/j^), solid carbon (Mq) and SiC (Mgic) are also shown. All the masses are expressed in solar units. 
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2.0 
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8.19d-06 
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Figure 10. Left: the variaton of t he surface cont ent of oxygen and carbon in two models of mass 6Mq in which the mass loss was modelled 
according to the prescriptions by iBlo ckei] {l99^ (solid, black line) and Vassiliadis & Wood (1993) (dotted, red). In this latter case only 
a little fraction of the stellar mass is lost by the time that the surface oxygen almost disappears, thus the star enters the C-star regime, 
and the silicates produ ction stops . The right panel (to be compared to the left panel of Fig. ??, that shows the corresponding evolution 
calculated with the lBlockerl lll995l'l recipe for mass loss) shows that the silicates production is extremely poor, and the star produces iron 
for the almost totality of its evolution. Note that the small carbon abundances inhibit the formation of C-rich dust particles, even in 
conditions of C/O > 1. 



relative to pyroxene and quartz (see Table 1 for the values 
adopted in the reference model). The resulting dust masses 
are shown in Fig. [TT] where dashed (solid) lines indicate the 
mass of pyroxene (quartz) expected when the minimum and 
maximum sticking coefficient is adopted. For comparison, we 
also show the total mass of silicates expected in the refer- 
ence model (filled squares) and mass of olivine (long-dashed 
line) which is not affected by the uncertainty on the stick- 
ing coefficient and is the most abundant species among sil- 
icates. Thus, although the pyroxene and quartz masses can 
vary by a factor ~ 20, their variation has a smaller effect 
on the total mass of dust: this is represented by the shaded 



region and shows that when the lowest limits for the pyrox- 
ene and quartz sticking coefficients are used, their contribu- 
tion is completely negligible compared to the olivine, caus- 
ing a ~ 30% decrease in the total mass of silicates. On the 
other hand, with the highest sticking coefficient, pyroxene 
becomes the dominant species for all stellar models, and the 
silicates mass produced increases by a factor ~ 3. At present, 
this quantifies the degree of uncertainty due to choice of the 
sticking coefficients on the total dust mass formed. 
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Figure 11. The filled squares indicate the total mass of silicates 
estimated by the reference model as a function of the initial stellar 
mass with the shaded region illustrating the uncertainty induced 
by variations of pyroxene and quartz sticking coefficients. The 
dashed (blue) and solid (red) lines indicate, respectively, the low- 
est/highest masses of pyroxene and quartz produced when the 
corresponding sticking coefficients are decreased/increased by a 
factor 5 with respect to the reference values shown in Table 1. 
The long-dashed line indicates the mass of the olivine— type dust 
produced in the reference model. 

6 DISCUSSION 

In Fig. [T2] we summarize the results of the present study 
of dust production by AGBs and SAGBs stars. We show 
separately the silicates (left panel) and carbon (right) dust 
produced by stars as a function of their initial mass. The 
values obtained with the reference model for the AGB evo- 
lution described in Sect.[2]are indicated as filled squares. The 
shaded areas quantify the uncertainties introduced in dust 
mass estimates by the poorly understood physical processes 
discussed in the previous sections. These are limited to stel- 
lar masses M 6Mq because a detailed investigation of the 
uncertainties affecting dust production by SAGBs will be 
presented in a future paper. For comparison, the tr i angles 
represent the results obtained by iFerrarotti fc Gail ()2006l ) 
using synthetic models of AGBs with initial metallicity of 
Z = 0.001. 

The extension of the shaded region in the left panel of 
Fig. [12] confirms that the mass of silicates produced by stars 
with M ^ 3A/0 is extremely uncertain, being very sensi- 
tive to the adopted molecular opacities whose treatment de- 
termines whether HBB can be activated. In our standard 
case, models within this range of masses are precited to pro- 
duce carbon-type dust, with a poor production of silicates: 
this explains the almost null lower limit in the left panel of 
Fig. 1121 On the other hand, use of c lassic opacities, coupled 
with the IVassiliadis fc Wood! (|l993l ) mass loss prescription, 
allows HBB conditions to be reached even for M ~ 2Mq, 
thus preventing the envelope of the star to become carbon- 



rich, and favouring the production of silicate-type dust. 
These assumptions, coupled with the highest choice of the 
sticking coefficients for the formation of pyroxene and quartz 
grains (see Fig. Ill[) . determine the upper envelope of the 
shaded area in the low-mass regime of Fig. [12] (left panel). 

For stars experiencing HBB, the uncertainties on sili- 
cates production are smaller. The lower limit on the shaded 
regions is set by the mass loss modelling, whereas the up- 
per limit is due to the choice of the sticking coefficients. In 
agreement with the discussion of the previous sections, we 
find that the most uncertain models are those experiencing 
strong HBB (M ~ &Mq). 

The right panel of Fig. [12] shows that carbon dust pro- 
duction is also uncertain in the low-mass regime (M^ ?,Mq), 
due to the combined effects of mass loss and opacities 
for carbon-rich mixtures, that are essential in determining 
whether HBB conditions are reached. In our reference model 
no carbon dust is expe cted for M > S.SMp, where as in mod- 
els calculated with the lVassiliadis fc WoodI (|l993l ) mass loss 
treatment some production is expected, because the enve- 
lope is consumed so slowly that eventually TDU can lead 
to the C/0> 1 condition. This determines the extent of the 
shaded region associated to stars of 4Af0 in the right panel 
of Fig. [H 

A comparison between the data marked with squares 
and triangles in Fig. [12] s hows that the dust yiel ds based 
on synthetic AGB models ijFerrarotti fc Galll2006l ) can not 
be reconciled with our results, even taking into account 
the uncertain ties on the physica l proce sses discussed so far. 
According to iFerrarotti fc Gail (l2006l ) , the dominant dust 
species produced by stars with Z — 0.001 and initial stellar 
masses IMq ^ M ^ 7Mq is carbon-dust; the production of 
silicates is limited to stars with initial masses ^ AMq but 
amounts to a few percent of the total dust mass formed. 

These differences can be entirely ascribed to the 
adopted description of the convective transport in the in- 
nermost regions of t h e conv ective envelope. As shown by 
I Ventura fc D'Antonal (|2005al ). use of the FST model leads 
to stronger HBB, at odds with the MLT treatment, that 
favours a modest HBB only in the most massive stars. In 
the latter case, independently of the mass, TDU eventually 
dominates and the carbon-st ar stage is reached. It i s not 
surprising that the results by iFerrarotti fc Gaill (|2006i ) pre- 
dict the formation of C-type dust; in our case this is possible 
only for models not experiencing HBB, i.e. for MsC 3Mq. 

In the low-mass regime the results are in qual- 
itative agreemen t , the carbon-dust mass presented by 
iFerrarotti fc Gaill li2006l ) being lar ger; this is partly due to 
the different choices for the low-T molecular opacities and 
the mass loss description, and also the efficiency of the TDU 
extension, which in our case is found by neglecting any over- 
shoot from the base of the envelope, and is thus to be con- 
sidered as a conservative estimate. 

6.1 Implications for cosmic dust enrichment 

It is interesting to investigate the implications of the dust 
yields for AGB and S AGB stars within the context of cosmic 
dust enrichment. 

Following the analysis done by Valiante et al. (2009), 
we can estimate the time evolution of the dust mass (and 
of individual dust species) produced by a stellar population 
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Figure 12. The mass of silicates (left) and carbon dust (right) formed around AGB and SAGB models of a given initial mass. Data points 
represented by filled squares indicate the results of our reference AGB and SAGB models. The shaded areas in both panels i llustrate the 
overall uncertainties associated to these predictions (see text). For comparison, the dust yields predicted bv lFerrarotti &: Gail (2006.) are 
also shown (filled triangles). Times on the upper axis indicate for each initial mass the duration of the evolutionary phases preceeding 
the beginning of AGB. 



as, 

Mdust(t) = / dt' mdnst{m)(j]{m)SFR{t' - Tm)dm, 

JO Jm,{t') 

(20) 

where mdust(?Ti) are the stellar mass dependent dust yields, 
4>{m) is the stellar initial mass function (IMF), rriup is the 
upper limit of the stellar mass range, Tm is the lifetime of 
a star with mass m, SFR is the star formation rate and 
m^it') is the mass of a star with lifetime Tm, = t' . We 
adopt a Larson IMF (Larson 1998) which follows a Salpeter- 
like power-law at the high-mass end but flattens below a 
characteristic stellar mass, rrich ~ 0.35Mq, 

(?i(m) ocm-^^+'^e"'"'^''/™, (21) 

where a = 1.35 and the IMF is normalized between minf = 
O.IMq and mup = WOMq. For stars with masses in the 
range IMq < m ^ 8Mq, we take the AGBs and SAGBs 
dust yields predicted by the reference model discussed in 
Sec. U At larger stellar masses, dust can be synthesized in 
the ejecta of core-collapse SNe and dust yields as a function 
of the progenitor mass and metallicities have been published 
in the literature (Kozasa, Hasegawa & Nomoto 1991; Todini 
& Ferrara 2001; Nozawa et al. 2003; Bianchi & Schneider 
2007). Here we neglect this contribution to dust production 
and we only consider the AGBs and SAGBs dust yields. 

The time evolution of the total dust mass, described 
by eq. 1201 depends also on the stellar lifetimes and on 
the adopted star formation history. For low and interme- 
diate mass stars, lifetimes are computed directly from the 
ATON evolutionary models and range between 72.4 Myr 
and 1.82 Gyr for AGBs and between 41.7 and 61.7 Myr for 
SAGBs. 

Fig. [13] shows the evolution of the total mass of dust 
(solid line) , silicates (dashed line) and carbon grains (dotted 
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time [Myr] 

Figure 13. Time evolution of the mass of dust produced by 
AGBs/SAGBs normalized to the total mass of stars (solid line). 
All stars are assumed to form in a single burst at time = 
with a metallicity of Z = 5 X 10^^ Zq and a Larson IMF with 
^ch = O.35M0. The dashed and dotted lines show the separate 
contributions of silicates and of carbon dust, respectively. 

line). All the stars are assumed to form in a single burst at 
t — with initial metallicity of Z = 5 x W^^Zq and the 
mass of dust is normalized to the total stellar mass formed. 

Since massive AGBs and SAGBs stars are the most 
efficient silicates producers, about 50 Myr after the burst 
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silicates dominate the dust mass evolution. When stars of 
M ~ 3M0 have ended their evolution, about 300 Myr af- 
ter the burst, the production of silicates ceases and restarts 
only when stars with M ^ 2Mq reach their AGB phase 
(~ 1.7 Gyr after the burst). Carbon dust enrichment is de- 
layed to longer timescales, being produced only by stars with 
meisses < 3M0; this represents the dominant dust compo- 
nent after ~ 470 Myr of the burst. 

Thus, if AGBs and SAGBs stars are the dominant star- 
dust sources, the ISM in galaxies at z > 10 and in young 
starbursts with ages < [450 — 500] Myr at any redshift will 
be predominantly enriched by silicates and no carbon dust 
features are expected to be observable. It is clear that these 
findings may partially depend on the adopted initial metal- 
licity of the stars, which affect both the core mass and the 
surface chemistry of the stars. In a future study, we will in- 
vestigate the formation of dust using a grid of AGBs/SAGBs 
exploring the full range of metallcities Q ^ Z ^ Zq . 

Presence of carbon dip in the (rest-frame) UV spec- 
trum of galaxies at z~ 1 - 3 (|No11 et al.ll2007l . 120091 ') suggests 
that these host evolved (age exceeding a few hundred mil- 
lion years) stellar populations. The lack of carbon dip in the 
spectra of AGN and high redshift quasars, accompanied by 
the presence of the silicate feature at 10/im, suggests that 
these active nuclei are hosted in young stellar populations 



('Maiolin o et~alll200ll : iLutz et al.ll2008l : iMaiolino et ai]|2004l : 
[Gallerani et al.ll2010l ). 



7 CONCLUSIONS 

We investigate the production of dust using physical models 
of AGB and SAGB stars with mass in the range IMq ^ 
M ^ 8M0 and initial metallicities Z = 0.001. 

The type of dust formed depends on the HBB ignition at 
the bottom of the external mantle: HBB is accompanied by 
the depletion of surface carbon, which, in turn, prevents the 
production of carbon dust, in favour of silicates. The thresh- 
old mass above which only silicates are formed is found to 
be M = 3M0. Below this limit, the evolution in the sur- 
face chemistry is dominated by TDU with the consequent 
formation of carbon dust. 

The distribution with stellar mass of the amount of dust 
formed shows a minimum corresponding to the [3.5 — 5]M0 
mass range due to the combined effects of the soft HBB 
experienced by these stars and inefficient mass loss. The 
strength of HBB increases with stellar mass leading to al- 
most complete destruction of surface oxygen for stars with 
M ~ 6M0; this, in turns, limits the formation of silicates 
and leaves solid iron as the dominant dust species. 

SAGB models are expected to produce great amounts 
of silicates. They experience large mass loss rates, ejecting 
their external envelope before significant surface oxygen de- 
pletion. 

In contrast with previous studies, we show that use of 
the FST treatment, as in the present investigation, leads to 
an efficient HBB ignition in all models with masses M > 
3M0, that are thus predicted to produce silicates. This is at 
odds with previous results published in the literature coming 
from synthetic AGB models based on the MLT treatment, 
which predicts C-rich dust to be the dominant species at all 
masses. 



As a consequence of the strongly mass-dependent dust 
composition found by our physical AGB models based on 
FST treatment of turbulent convection, whenever the ISM 
dust enrichment is dominated by AGBs and SAGBs at low 
metallicities, we predict an early phase (~ 300 Myr) of sili- 
cate production followed by the formation of C-rich dust. 
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